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Abstract 

We calculate the shadowing of sea quarks and gluons and show that the shad- 
owing of gluons is not simply given by the sea quark shadowing, especially at small 
x. The calculations are done in the lab frame approach by using the generalized 
vector meson dominance model. Here the virtual photon turns into a hadronic 
fluctuation long before the nucleus. The subsequent coherent interaction with 
more than one nucleon in the nucleus leads to the depletion a{^* A) < Aa(^*N) 
known as shadowing. A comparison of the shadowing of quarks to E665 data for 
Ca and 207 Pb shows good agreement. 
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1. Introduction 

When calculating perturbative QCD cross sections in nucleus nucleus 
collisions one has to take care of an additional effect not appearing on the 
pure nucleon nucleon level: nuclear shadowing. As a result of this depletion 
of the nuclear parton densities at small x one finds a strong suppression of 
e.g. charmonium states and minijets 0. Also for prompt photons smaller 
multiplicities result due to the smaller number density of partons in the 
relevant region of the momentum fraction variable. 

In either model (lab frame vector meson dominance type models or infinite 
momentum frame parton fusion models) the shadowing of gluons is expected 
to be much stronger than the shadowing of sea quarks, even at small x. This 
seems to contradict the first naive expectation in terms of the QCD improved 
parton model that explains the scaling violation of the structure functions via 
the DGLAP splitting functions which treat sea quarks similar to gluons (in 
the sense that the sea quarks are produced by the gluons when Q 2 increases) 
in a region where essentially no momentum is carried by valence quarks, 
i.e. in the small x region. 

In the following we will focus on the lab frame interpretation which explains 
the shadowing phenomenon by use of the generalized vector meson dominance 
model (GVMD). 

2. Lab frame description of shadowing 
In the lab frame description one essentially makes use of the hadronic 
structure of the virtual photon, manifesting itself in a field theoretic approach 
in two different time orderings assuring gauge- and Lorentz invariance to the 
amplitudes At small enough x (x <C 0.1) the hand-bag graph contribution 
gets small. The interaction then proceeds via the VMD graph where the 
virtual photon fluctuates into a qq pair within the "coherence time" A « 
2v/ {fJ? + Q 2 ) ~ l/(2mx), where m is the nucleon mass, /i the mass of the pair 
and v the lepton energy loss. At small Q 2 the interaction with a free nucleon 
then proceeds via the p,u, and mesons. At larger Q 2 the qq continuum 
has to be taken into account to describe the F^(x,Q 2 ) data. Including sea 
quarks xfjy(x,Q 2 ) only, the nucleon structure function can (for small x) be 
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written as 



F 2 N (x,Q 2 )=^xf N (x,Q 2 ) = (1) 

With the transverse size of the qq pair frozen during the scattering, due to 
Lorentz contraction, the photon-nucleon cross section can be factorized: 

a( 1 *N) = £dz J d 2 r\^(z,r)\ 2 a qq - N (r) (2) 

with the Sudakov variable z. The photon wave function \ip(z, r)\ 2 can be 
interpreted as the propability of finding a qq state with transverse separation 
r and a momentum fraction z and (1 — z) with respect to the virtual photon. 
The wave function can be expressed as || 

\^(z,r)\ 2 = ^£e^(*)« a *i(«f (3) 
z " i=i 

where P(z) = z 2 + {l — z) 2 is the splitting function of the gauge boson into the 
qq pair, a 2 = Q 2 z(l — z) and ii'i is the modified Bessel function. The most 
important contributions come from the region where ar is small. Therefore 
pairs with a small transverse size are favored. 

The cross section for the strong interaction of the pair with the nucleon, 
which happens via gluon exchange, is 

P q qN = ^r 2 a s {Q' 2 )x'g N {x\ Q' 2 ) (4) 

with x' = x/a 2 r 2 and Q' 2 = 4/r 2 . The strong scaling violation of x'g^{x' ', Q' 2 ) 
at small r (large Q' 2 ) in turn compensates for the smallness of the pair (for 
further details see |J and references therein). Therefore the exact expression 
for the sea quark density is 



xf N (x,Q 2 ) = -^-rx [ f ^--a q q N (x',r 2 ) (5) 
A-k 6 Jx x' 2 J4/Q 2 r 4 



The main problem at this point arises from the interplay between hard and 
soft physics. In some region A/Q 2 < r 2 < 4/Qq perturbation theory should be 
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valid. For r 2 > 4/Qq, i.e. for small Q 2 , the non-perturbative part dominates. 
To use the presently available parametrizations of the free nucleon parton 
densities we choose a lower cut-off of Ql = 0.4 GeV 2 . It is obvious that 
we thereby completely neglect the non-perturbative input below this scale. 
However, in 0] it was shown that for Q 2 < 10 _1 GeV 2 a saturation in the 
shadowing ratio sets in. As a result, the uncertainty by choosing a cut-off 
at Qq = 0.4 GeV 2 is on the few percent level. Also, when comparing to 
experimental results, we find good agreement with the data. 
In the case of the gluon distribution of the free nucleon one considers the 
scattering of a virtual colorless probe (e.g. the virtual higgs boson) which 
proceeds via the production of a gluon pair which then strongly interacts 
with the nucleon. With the momentum cut-off the gluon density can be 
written as 



The main difference between the scattering of the qq pair and the gg pair 
stems from the much larger cross section cr ggN = 9/4 a q q N which in turn leads 
to a much stronger gluon shadowing. 

In the case of deep inelastic scattering off nuclei the same relations hold 
but with the respective cross sections for the scattering of the hadronic and 
gluonic fluctuations off the nucleus. For the sea quark distribution in the 
nucleus A one has 



As stated above, it is the long distance phenomenon of the hadronic fluctu- 
ation which causes the shadowing effect. For the small x, i.e. large A, the 
fluctuation interacts with more than one nucleon inside the nucleus. As a 
result of this coherent interaction one finds that a^A) < Aa(j*N). The 







and for the gluon the resulting equation is given by 
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specific feature that is responsible for the shadowing effect is the fact that 
(TfiA 7^ Aa hN . Within Glauber theory in the eikonal limit one finds [^] 

a hA = 2j d 2 b (l - e -™^(6)/2) (g) 
with the nuclear thickness function 

T A (b) = I dzp A (b, z) = -^e-" 2 ^ (10) 

The integration can be done exactly for a Gaussian parametrization of T A (b). 
Here we will include only the double scattering contribution which was shown 
[§] to strongly dominate the overall shadowing ratio Rq = 9n(x, Q 2 )/g A ] the 
triple and higher scattering terms give a contribution that is only on the 
percent level and is therefore neglected here. We use @, |8j 

a hA = Aa hN \l-Al exp j + . . . 

with a = 1.1/m. 

3. Results 

Based on the approximations above, namely double scattering contribu- 
tion only, a lower momentum cut-off at Ql = 0.4 GeV 2 , due to the lack 
of information on the non-perturbative input, and by employing the Gliick, 
Reya, Vogt parametrization we derived the shadowing ratios for quarks 
and gluons at the typical semi-hard scale Q = 2 GeV for 207 Pb and 40 C a (see 
figures [j] and ^j. One clearly sees that gluons are much stronger shadowed 
than sea quarks at small x which is due to the effects outlined above. This 
is a very important feature since at the soon available collider energies one 
will particularly test the small x region (x = pr/y/s at midrapidity) which 
becomes increasingly important due to the strong gluonic component inside 
the free nucleon and therefore in the heavy nuclei. 

This feature has a large impact on transverse energy production and on the 
initial temperature and entropy produced in heavy ion collisions. It should 
therefore come as no surprise when the final particle multiplicities are much 
lower than expected from calculations using e.g. the quark shadowing ratio 



(11) 
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Figure 1: R F2 vs. R G at Q 2 = 4 GeV 2 for 207 Pb. 



6 



1.2 




Figure 2: R F2 vs. R G at Q 2 = 4 GeV 2 for 40 Ca. 
Rp 2 {x) shown in figure || for gluons also. To check our results, we compared 



the results for Rp 2 (x) to Fermilab-E665 data [1C] and find a good agreement. 
Since s ~ Q 2 /x, for deep inelastic events, one is restricted to a certain Q 2 - 
range for the different x- values in the experiment. To account for this feature 
we compared to the data within the range (Q 2 ) ~ 2.42GeV 2 - 4.45GeV 2 be- 
cause our calculation was done at fixed Q 2 = 4 GeV 2 . One also sees that the 
agreement with the data is only good at small enough x, i.e. at large enough 
coherence lengths A ~ 1/Ak z , where Ak z = kj — k^ is the phase shift between 
the photon and the hadronic fluctuation. To derive the shadowing ratio for 
the complete structure function F 2 one could in principle improve the large x 
behavior by inserting a splitting function P gq = z 2 + (1 — z) 2 as suggested in 
U where z = x/x' . Also one has to include the valence quark distribution as 
F^(x, Q 2 ) = J2i e 2 [x/n^x, Q 2 ) + xVi(x, Q 2 )]. But we here neglect this point 
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Figure 3: R F2 = F^/AFJ? as parametrized by Eskola in @] at Q 2 = 4 GeV 2 . 

since we were mainly interested in the region of small x. 

4. Conclusions 

We calculated the shadowing of sea quarks and gluons at small values 
of x in the generalized vector meson dominance model (GVMD) approach. 
Based on the Glauber-Gribov multiple scattering picture for a qq pair with 
a coherence length 1/Ak z = A ~ l/2mx, we showed that one should expect 
significant differences between Rp 2 and Rc at small x. This feature, often 
neglected in the past, has severe consequences on quantities such as char- 
monium production, minijets and therefore on the production of transverse 
energy and entropy When applying the much stronger gluon shadowing 
one should expect significantly smaller multiplicities for processes involv- 
ing gluons in the initial state compared to calculations that accounted for 
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all the shadowing effects by simply employing the quark shadowing ratio 
Rp 2 = F^/AF^ for all parton species. 
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